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1 Introduction 


Over the past years significant progress has been made in the determination of the strong 
coupling constant a s . Next-to-leading order (NLO) theoretical predictions are generally 
available. For some inclusive quantities also the next-to-next-to-leading (NNLO) orders have 
been calculated and estimates of the next higher terms exist. In addition, resummations 
of leading-log (LL) and next-to-leading-logarithmic (NLL) corrections to all orders have 
been performed in some cases. Power law corrections are controllable in the framework of 
the operator-product-expansion (OPE) or the resummation of renormalon chains. On the 
experimental side information exists from many different reactions over a large range of both 
space-like and time-like momentum transfers. Reactions include neutrino- and lepton-nucleon 
scattering, proton-(anti)proton collisions, e + e _ -annihilation and decays of bound states of 
heavy quarks. Observables are total cross section measurements, sum rules, scaling violations, 
branching ratios, global event shape variables and production rates of hadron jets, i.e. bundles 
of particles close by in phase space originating from isolated hard partons. 


1.1 Theoretical Predictions 


The QCD prediction for a cross section a(Q) at an energy scale Q can be expressed as sum of 
perturbative terms d per t, which vary logarithmically with energy, and non-perturbative power 
law corrections 5 np . The perturbative part is of the form 

fipert — CX s (fJj)A + Oi 2 s (n) ^ B + A — In + ... 

Here /i is an arbitrary renormalization scale used in the calculation and /3q the leading order 
coefficient of the QCD beta-function. For n ~ Q the perturbative expansion is well behaved, 
allowing for a measurement of the strong coupling constant. For |ln(p 2 /Q 2 )| 1 the 

effective expansion parameter becomes aq(/i) ln(/i 2 /Q 2 ), which leads to an unstable theoretical 
prediction. The full theory does not depend on p, but the truncated perturbative expansion 
does. Variation of /j thus allows to probe the sensitivity of the theory to uncalculated higher 
orders and thereby to assess theoretical uncertainties. 

In the following the nominal scale /i of an a. 5 -measurement is identified with the physical 
scale of the respective process. For data covering a range of scales the geometric mean of high 
and low end is taken, unless a central scale is specified explicitly by the authors. Combined 
results from different measurements are quoted at the geometric mean of the individual scales. 
The geometric mean for the scales is motivated by the fact, that (in leading order) a s is inversely 
proportional to the logarithm of the scale. 

The beta function describes the renormalization scale dependence of the strong coupling, 
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with fa = 11 — 2n//3 and Pi = 51 —19u//3 for QCD based on an SU(3)-colour gauge symmetry. 
The quantity nf denotes the number of active quark flavours. Measurements of a s obtained 
at different energy scales can be compared, either by evolving backwards to the point A where 
a s diverges, or by evolving to a common reference energy, which in recent years has become 
the Z mass. In the evolution [[lj a s is continous at flavour thresholds, which, within the still 
comparatively large uncertainties of the quark masses, can be taken at the MS-running mass of 
the respective quark flavour. As an example consider the evolution of an a s -measurement with 




71/ = 3 from the scale p = m T to p, = Mz . It proceeds by first going back to the charm-quark 
mass m c = 1.3 ± 0.2 GeV using the evolution equation for three flavours, then going up to 
rrif, = 4.3 ± 0.2 GeV with rif — 4 and finally evolving all the way up to p = Mz with 7 i/ = 5 
active flavours. Many slightly different schemes are in use to perform the evolution numerically. 
Despite being formally equivalent to NLO the resulting differences in the value for a s (Mz ) can 
be as large as 5a s (Mz ) = 0.001. To be consistent, for the numbers presented in this paper the 
treatment of flavour thresholds described above has been applied together with Runge-Kutta 
integration of the NNLO-beta function. 

1.2 Combination of Individual Results 

Many measurements of the strong coupling constant are dominated by theoretical uncertainties. 
In order to be quantitative, it will be assumed that they represent 68 % confidence intervals. 
Their non-statistical nature suggests that they should be interpreted in the Bayesian sense. Note 
that doing this also for the experimental uncertainties would justify to combine experimental 
and theoretical errors in quadrature. 

Theoretical errors for different ay-measurements are correlated via the underlying theory. 
Such correlations could be quantified by the derivatives of the predictions with respect to the 
various sources of uncertainty, e.g. unknown higher order coefficients. Instead, in many cases 
only total errors are available, leading to a situation where measurements are known to be 
correlated with very little information about the actual size of those correlations. To deal with 
this sitution, the following averaging procedure f 2 | has been employed. 

For a set of measurements Xi ± Uj the average is calculated as x = Y^ x i w i, with weights 
inversely proportional to the squares of the errors a*. This is a robust estimate of a common 
mean, which is also optimal if the single measurements are uncorrelated. To determine the 
error of the average and its y 2 the full covariance matrix is needed. Using only the diagonal 
terms, both are underestimated in the presence of positive correlations, i.e. a y 2 smaller than 
its expectation value is indicative of such terms. In this case a realistic error for the average 
can be obtained even if the size of the correlations is not known a priori, by introducing off- 
diagonal terms poycr/ in the covariance matrix and adjusting the effective correlation coefficient 
p such that y 2 attains its expectation value. Averages with errors determined according to 
this method will be referred to as “conservative averages” in the following. Measurements with 
asymmetric errors are shifted to the center of the error band and the error symmetrized prior 
to the averaging procedure. 


2 Measurements of a s 

This section contains an overview over the most important types of ay-measurements, presenting 
results from a large variety of reactions and observables. 

2.1 Inclusive Measurements 

Inclusive measurements of ay depend only on one energy scale characterizing the process and are 
best understood theoretically. The perturbative prediction is known to (P(a^), non-perturbative 
effects can be treated in the framework of the Operator Product Expansion ||(OPE). 


The Ratios R 7 and Ri 

The quantities R 1 and Ri are defined by the ratios of the hadronic to the leptonic branching 
fractions of a virtual photon or the Z-boson. In both cases the hadronic system is formed from 
the electroweak (EW) coupling of a vector boson (Photon or Z) to a primary quark-antiquark 
pair. The sensitivity to the strong coupling comes from gluon radiation off the primary quarks. 
This radiation opens up new final states for the hadronic system, which increase the hadronic 
width with respect to the purely electroweak expectation. The QCD correction is dominated by 
the perturbative terms. The leading non-perturbative effects scale with 0(l/Q 4 ) and thus are 
negligible for and Ri. Mass effects, in particular from the large top-bottom mass splitting, 
have been calculated Q to 0(a 2 ). 

The first analyses || of were still based on an erroneous third order correction ||. A 
recent reanalysis [j?] using the correct third order coefficient [3'] and data from centre-of-mass 
energies between 5 GeV and 65 GeV yields a s (/x) = 0.175±0.023 for a central scale n — 18 GeV, 
which corresponds to a s (M z ) = 0.128±§;§i§. 

A simple parametrization [;|] for R{ is available to third order in a s (M z ) . Using the 
combined result from all four LEP-experiments [10| Ri = 20.778 ± 0.029, a top-quark mass 
m t = 180 GeV/c 2 and a Higgs-mass M H = 300 GeV/c 2 one finds a s (M z ) = 0.124 ±0.005. The 
statistical error is 5a s (M z ) stat = 0.004, twice as large as the QCD theoretical uncertainty. 



Figure 1: Combined measurement 
of a s and the Higgs mass M H from 
QCD radiative corrections to the 
standard model. 


In addition to Ri also the hadronic peak cross section of the Z and its width are sensitive 
to the strong coupling. Combining this information with all available electroweak data from 
LEP, SLC, collider measurements and Deep Inelastic Scattering in a global standard model fit, 
allows the simultaneous determination of a s (M z ) and the Higgs mass. The result is shown in 
Fig. [I]. One finds [|TT| a Higgs mass M H = 149 GeV/c 2 and a s (M z ) = 0.1202 ± 0.0033. The 
standard model constraint thus not only gives a very precise measurement of a s (M z ), but also 
one perfectly consistent with other precision measurements. 









Measurement of a s from R r 


An ci! s -measurement can also be obtained from R T = D hadr / B e , the ratio of the hadronic to the 
electronic branching ratio of the tau lepton. Here QCD radiative corrections affect the hadronic 
final state from a W-decay. In contrast to a s -determinations from or Ri the mass of the 
hadronic system is not fixed in tau decays. This makes the quantity Rr double inclusive, i.e. 
integrated over all hadronic final states at a given mass and integrated over all masses between 
M n and m T . 

The low energy scale requires a good understanding of the non-perturbative contributions, 
which in the framework of the OPE are proportional to vacuum expectation values 
(condensates) of the QCD fields. These condensates can be determined from independent 

B 


phenomenological analyses 


or together with a s from the higher moments of the mass 


spectrum of hadronic r-decays. B It turns out that the non-perturbative corrections to Rr 
are surprisingly small, S np = —1.5 ± 0.4%, and |14[ that perturbative QCD is applicable down 
to mass scales below 1 GeV. As a consequence an a s measurement from R T is potentially very 
accurate. 

Assuming the validity of the completeness relation for the tau branching ratios into hadrons, 
electrons and muons, B^ a dr + B e + B^ = 1, and lepton universality, the ratio Rr can be expressed 
as function of B e alone. Note that the larger mass of the muon leads to B p /B e = 0.9726. 
The branching ratio B e can be determined by direct measurements or, again assuming lepton 
universality, by comparing mass and lifetime of the tau lepton and the muon. From recent 
measurements JT5| one obtains Rr = 3.642 ± 0.010. A value for a s {m T ) with a rather 
conservative error based on this number is quoted 0 as a s [m T ) = 0.33 ± 0.03. This error 
is almost entirely due to uncertainties in the perturbative prediction S pert , taken to be half 
the difference between Le Diberder-Pich resummation 11161 and resummation of renormalon 


chains. [[17]] Other studies [[L8| suggest that it could be significantly smaller. 

A combined analysis of Rr and the leading moments of the hadronic mass spectrum has 
been performed by the ALEPH [0] and CLEO [|20| collaborations. A three-sigma discrepancy 
between the two was due to the old 1994 PDG-value [|0 for B e used in the CLEO analysis. The 
moments extracted from the hadronic mass spectrum are in good agreement. With an updated 
leptonic branching ratio the CLEO result a s (m T ) = 0.339 ±0.024 becomes compatible with 


the ALEPH number a s (m T ) = 0.353±0.022. The conservative average of both ^^-measurements 

HI 


is a q im T 


= 0.347 ± 0.022 ± 0.03, where the second error 


has been added to account for 


the fact that both analyses are based on Le Diberder-Pich resummation. Evolving up to the 
Z-mass and symmetrizing the larger one of the slightly asymmetric errors, one finally obtains 
a s {M z ) =0.122 ±0.005. 


a s from the Gross-Llewellyn-Smith Sum Rule 


Determinations of the strong coupling constant based on sum rules are fully inclusive 
measurements at a very low Q 2 -scale. In the Quark-Parton-Model (QPM) the Gross-Llewellyn- 
Smith sum rule (GLSR) counts the number of valence quarks in the nucleon. The perturbative 
QCD-correction |23|] is known to O(al). A measurement of the strong coupling constant taking 
into account also non-perturbative (“higher twist”) terms [^4[ has been performed by the 
NuTeV-Collaboration [^] at a scale ji = 1.73 GeV, based on the old CCFR data.||26 


result ct s (/i) = 0.260to;o 4 g corresponds to a s (M z ) = O.llOl^ddg. As the input data have been 
re-calibrated recently, this number can be expected to be updated, too. 

































































The Bjorken Sum Rule 


Also for the Bjorken sum rule (BjSR) the perturbative QCD-correction to the QPM [p? 


is 


known to O(al) together with an estimate of the size of the non-perturbative effects. |28fl The 
BjSR is defined for charged-current neutrino-nucleon scattering or polarized lepton nucleon 
scattering. Measurements of the strong coupling constant so far only exist based on the spin- 
dependent structure functions of the nucleon. A first analysis [|29| had comparable experimental 
and theoretical errors, a recent update [^Oj finds a s (/i) = 0.320^o;o53(exp) ± 0.016(theo) for a 
scale p = 1.73 GeV. Evolved to the Z-mass one obtains a s {M z ) = O.llStow with an error 
dominated by the experimental uncertainties. 

The small size Sa s (/j,) = 0.016 of the theoretical uncertainties results from the stability of 
the analysis with respect to various ways of estimating missing higher order terms by means 
of Pade approximants (PA). The PA [N/M] of a function is the ratio of two polynomials of 
order N and M, which to order N + M has the same Taylor-expansion as the original function. 
Pade approximants thus offer a systematic way to guess how a perturbative expansion resums, 
by rewriting a perturbative series as a ratio of two polynomials. Compared to the original 
expression, the PA introduces poles in the coupling plane in a similar fashion as expected 
from renormalons. This may explain why PAs seem to be able to approximately resum the 
perturbative series, a Ending corroborated e.g. by the observation |3l| that measurements of 
the strong coupling from global event shape variables become much more consistent when the 
second order prediction is replaced by its [1/1]-PA. 


Heavy Flavour Thresholds 

Another potentially very precise a s -measurement is derived from the threshold behaviour of 
the 66-production cross section in e + e _ -annihilation. QCD sum rules for this process are 
dominated by the non-relativistic threshold behaviour and allow to extract simultaneously a s 
and the 6-quark mass. A determination of the strong coupling at the BLM-optimized [[£|] 
renormalization scale /i = 3 GeV yields a s (M z ) = 0.110 ± 0.001. Naively one would set the 
renormalization scale to the 6-mass, which corresponds to a shift [|34| Sa s (Mz ) = +0.008. 
Taking this as the theoretical error of the measurement one has a s (Mz ) = 0.110 ± 0.008 or 
aq(/i) = 0.217 /q;o3o for p = 3 GeV, respectively. 


2.2 Bound States of Heavy Quarks 

The strong coupling constant has been determined from ratios of decay rates, which are 
described by perturbative QCD, and level splittings between different radial excitations of 
the bound system, which probe the QCD-potential between quark and antiquark. 


Decays of Heavy Quarkonia 


A precise a s -measurement is based on the ratio of the hadronic to the leptonic width of a heavy 
quark-antiquark pair. At Born-level this ratio of annihilation rate into three gluons over the 
rate into a lepton pair is proportional to O(o? s /a 2 ern ). Taking into account relativistic corrections 
proportional to the average < v 2 /c 2 > of the quarks, the theoretical prediction can be written 
in the form 


P (qq—* hadrons) 
T(qq—>e + e~) 



The perturbative prediction R per t is known to NLO, D is a free parameter for the size of the 
non-perturbative relativistic corrections. Assuming D to be a universal constant, it has been 














extracted |35j together with a s in a combined analysis of T and J/T decays. Theoretical 
uncertainties in R per t were studied by introducing ad hoc NNLO terms, by varying the 
renormalization scale and by Pade-like rewriting terms of the type (1 + Ba s ) as 1/(1 — Ba s ). The 
result at a scale /i = 10 GeV is a s (p) = 0.167to;on, dominated by the theoretical uncertainties. 

The CLEO collaboration |36| has extracted a measurement of the strong coupling from the 
ratio T(T— ygg^/T^T^ggg), which to leading order is proportional to a em /a s . The result 
at the scale of the Upsilon mass p = 9.7 GeV is a s (p) = 0.163 ± 0.009(exp) ± O.OlO(theo). 
The experimental error is dominated by systematic uncertainties in the photon background 
of purely hadronic Upsilon decays and roughly equal to the renormalization scale error. The 
conservative average of both measurements taken at the scale p = 9.7 is ct s (p) = 0.166 ± 0.013, 
which evolves to a s (Mz) = 0.112 ± 0.006. 


Lattice Calculations 

Very precise determinations of a s were performed in the analysis of level splittings between the 
S- and the P-states of the T-system by means of lattice calculations. |I7| The current level of 
precision is a result of reduced lattice spacing errors, a better understanding of the conversion 
from the lattice to the MS coupling measured elsewhere and the introduction of fermion loops 
on the lattice. Calculations exist with rif — 0 and rif = 2 dynamic fermions, which give 
only marginally different results and thus allow a safe extrapolation to the physical case of 
rif = 3 light flavours. The most recent results [|38| are a s (M z ) = 0.118 ± 0.003 (NRQCD) and 
a s (M z ) = 0.116 ± 0.003 (FNAL). The conservative average is a s (Mz ) = 0.117 ± 0.003, which 
at the typical scale of the calculations p = 8.2 GeV corresponds to ct s (p) = 0.184 ± 0.008. 


2.3 Scaling Violations 

Scaling violations are observed in Deep Inelastic Scattering (DIS) processes in the structure 
functions of the nucleon, and in e + e _ -annihilation into hadrons in the fragmentation functions of 
the primary partons, i.e. in reactions with space- and time-like momentum transfers. In both 
cases perturbative QCD predicts a softening with increasing Q 2 , d In F/d In Q 2 cx a s (Q) as 


described by the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations. |39 


The softening of structure functions comes about because higher momentum transfers resolve 
more partons from vacuum fluctuations in the nucleon. For fragmentation functions the growing 
phase space permits additional gluon radiation and thus enhanced particle multiplicities in 
jets. The theory |4(J is known to next-to-leading order. In a determination of a s perturbative 
and non-perturbative effects can be disentangled through their energy dependence. Here DIS 
processes are favoured because non-perturbative effects decrease rapidly with 1 /Q 2 . In e + e~- 
annihilation they decrease only proportional to 1/Q, which renders measurements of a s from 
fragmentation functions less precise than the ones from DIS, despite the fact that they are 
performed at larger momentum transfers. 


Scaling Violations in Structure Functions 

Measurements of the strong coupling constant were performed in DIS with neutrino beams 

The most recent results for A^-L 
are collected in Tab. [I] and displayed in Fig. |j. All 
numbers are in good agreement, although there is a slight trend towards larger values for more 
recent measurements. The two most precise results are from a combined analysis of SLAC- 


and charged leptons on targets of heavy and light nuclei, 
from various experiments HI [4|. |a 


BCDMS measurements W1 


A^ = 263 ± 42 MeV, and a recent re-analysis of CCFR data [[73 
























Measurement 

H/ GeV 

Agi/MeV 

CDHS 

uFe 

4.5 

250 ± ff 0 

CHARM 

uCa 

3.9 

310 ± 157 

BEBC 

z/Ne 

4.0 

100 ± If 

BCDMS 

ytC 

8.4 

230 ± 63 

EMC 

pH 2 

4.7 

011 - 1 - 117 
^ 108 

SLAC-BCDMS 

HWD 2 

7.1 

263 ± 42 

CHDSW 

uFe 

6.7 

300 ± 100 

NMC 

pH 2 ,^D 2 

2.6 

306 ± li 

CCFR 

uFe 

8.2 

346 ± 58 

Average 

5.4 

287 ±31 


Table 1: Measurements of the QCD scale from DIS in chronological order. The errors are the purely 
experimental uncertainties. The average is the weighted average from the lower part of the table. 


which gives A^b = 346 ± 58 MeV. For the latter the logarithmic derivatives of the structure 

The value for moved 
based on the same data, mainly as a 


functions together with the NLO QCD-fit are shown in Fig. 
up by 136 MeV compared to a previous publication 


26 


consequence of an improved energy calibration for the detector. The weighted average based on 
the measurements published after 1990 yields A^-A = 287±31 MeV. With a common theoretical 


uncertainty [42| of 5a s (M z ) = 0.004 this corresponds to a s (M z ) = 0.116±0.005, or equivalently 
a s (/j,) = 0.200 ±0.016 for an average scale yt = 5.4 GeV. The total error from this NLO analysis 
is still dominated by theoretical uncertainties. First results from NNLO-analyses [45| indicate 
that the value of a s is quite stable when going to higher orders, i.e. the theoretical error can 
be expected to decrease. 
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Figure 2: Graphical representation of the 
results given in table [j]. The vertical line 
marks the average A^, the shaded area its 
uncertainty. 
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Figure 3: Logarithmic derivatives of the structure functions F 2 and xF 3 as 
measured by the CCFR collaboration. The solid line is the result of a NLO 
QCD-fit. 


Scaling Violations in Fragmentation Functions 


The full NLO theoretical framework has also been used in two determinations of the strong 


coupling constant from scaling violations in fragmentation functions. ||46| , |47[ Both are based 
on inclusive and uds, c, b and gluon-jet enriched data samples from hadronic Z decays, 
combined with measurements from lower centre-of-mass energies down to yfs = 22 GeV or 
y/s = 14 GeV, respectively. The strong coupling constant was determined together with 
parametrizations for the fragmentation functions of the different parton types and a power- 
law correction describing non-perturbative effects. The latter was found to be small for the 
energy range under consideration. The results of these model independent analyses are ph) 
a s (M z ) = 0.126 ± 0.009 and [£T7| a s (M z ) = 0.121 ± 0.012, where the error is dominated by the 
QCD scale uncertainty. The larger error from the second analysis reflects a larger range for 
scale variations. The conservative average of both results is a s (M z ) = 0.124 ±0.010, which for 
a central scale p — 36 GeV corresponds to a s (p) = 0.146 ± 0.014. 

Assuming that non-perturbative effects and differences in the fragmentation between quarks 
of different flavours and gluons are described correctly by Monte-Carlo models, one can also 
exploit the LUND matrix element model as the theoretical basis for a measurement of a s from 
scaling violations. A result |48| obtained by the DELPHI collaboration is consistent with the 
other determinations. Finally, it is worth noting that scale breaking effects in fragmentation 
functions can also be studied in DIS, using the momentum distribution of particles in the 
current jet.]|9| DIS thus offers the possibility to measure a s from scale breaking effects both in 
structure functions and in fragmentation functions. 


2.4 Results from Hadron Colliders 

Prompt photon production in hard parton-parton scattering is a Compton-like process of 
0(a s a em ). In the difference a(pp—^X) — a(pp^>'yX) the sea quark and gluon structure 
functions of the nucleon cancel, i.e. only the well known valence quark distribution is 


needed as external input for an a s -determination. The UA6 collaboration |5D] finds A^ = 


235 ± 106(exp)lg 46 (theo) MeV. Translated into a s at the typical scale p = 4 GeV of the 
measurement one finds a s (p ) = 0.206±o;o33, which corresponds to a s (M z ) = 0.112±o;[uo- 

Heavy quarks not present in the initial state can be produced by quark-antiquark 
























































annihilation or gluon-gluon fusion processes, which to leading order are of 0(a^). The theory is 
developed to NLO. Experimentally those reactions can be tagged by the decay characteristics of 


the heavy hadrons. The most precise measurement |5l| from an analysis of bb +jets production 
in pp collisions is a s (/i) = 0.138±o;oi9, determined at a scale p = 20 GeV. Evolved to the Z 
mass one finds a s (M z ) = O.lOQ/o/j®. 

The strong coupling constant has also been determined from QCD-radiative corrections to W 
production at hadron colliders, based on the cross section ratio R w = a(W+l jet)/cr(iy+0 jets). 
This ratio is sensitive to the product a s F(x ) of strong coupling constant and parton densities in 
the nucleon. The perturbative QCD-correction is known to NLO. Taking the nucleon structure 
functions from low energy data, a s enters both in the matrix element for jet-production in 
W decays and in the evolution of the structure functions to collider energies. A measurement 
a s (M w ) = 0.123 ± 0.025, or equivalently a s (M z ) = 0.121 ± 0.024, has been obtained by the 


UA2-collaboration [^] at a centre-of-mass energy a/s = 630 GeV, with an error dominated by 
experimental uncertainties. The ratio Rw has also been determined at the Tevatron [[53]], but 
it turns out to be rather insensitive to a s when parton densities from low energy measurements 
are employed in the analysis. The reason is that for larger values of a s the enhanced softening 
of the structure functions in the evolution almost compensates the increased jet production 
rates from the matrix element. In other words, the determination of the strong coupling from 
W-decays at the Tevatron requires parton densities measured at the same energy. 

Another interesting prospect is the determination of the strong coupling constant from the 
inclusive jet cross section at hadron colliders. Here an extraction of a s (-EV) with high precision 


seems feasible. [54 


2.5 Global Event Shape Variables 

To leading order the fraction of three-jet events in e + e~-annihilation is proportional to a s . 
Global event shape variables which are sensitive to the topology of multijet events exploit this 
fact for a measurement of the strong coupling constant. In order to be defined in perturbative 
QCD they must be constructed such, that they are insensitive to soft gluon radiation and 
that they remain unchanged if any of the final state particles splits into two collinear ones. An 
example is the variable thrust |55[ T = max r j(E p |p-n|)/(E P |p|), which measures the collimation 
of the momentum flow in an event. An ideal two-jet event without final state gluon radiation 
has T = 1. Gluon radiation decreases the value of T until one finds T = 1/2 for a perfectly 
isotropic event. 

Another example is y 3 , defined by means of a jet clustering algorithm where initially each 
particle is considered its own jet. Then those two jets which are closest in phase space are 
combined by adding their 4-momenta. Iterating the procedure, y 3 is defined as that distance 
where the event makes the transition from three to two jets. Common measures for the distance 


between two jets i and j are the Durham-metric ||56 
Jade-metric 


,d 


57 


Vij ~ 


ypj = 2 min (A/, EJ)(l — cos%)/s or the 


= 2(EiEj)(l — cos 6ij)/s. In both cases 8 l] is the opening angle between 
the jets and s denotes the total invariant mass of the hadronic system. 

The theoretical prediction for all global event shape variables is known to NLO, based on 
a numerical integration |58| of the ERT-matrix elements. |k]] For some variables also leading- 
logarithmic and next-to-leading logarithmic corrections, which dominate the cross section at 
high thrust, have been resummed to all orders, pm I 11 the latter case an improved theoretical 
prediction is obtained by combining the resummed prediction with the second order matrix 
elements, which is exact to O(a^) over the whole phase space and contains the dominant 
terms for two-jet like configurations to all orders. There is a certain freedom in performing the 
matching, |6lH which gives rise to differences at O(al). It thus can be employed as an alternative 






















































to the variation of the renormalization scale to probe the sensitivity of an a s measurement to 
unknown higher order perturbative corrections. 

The non-perturbative transition from partons to hadrons gives rise to power-law corrections, 
which in contrast to the case of inclusive variables cannot be handled in the framework of the 
OPE and thus usually are estimated by means of Monte Carlo models. |62| This dependence 
on phenomenological models introduces an additional “hadronization uncertainty” into an ad¬ 
measurement. Only recently analytic calculations |)3j became available, based on the study 
of the long-distance behaviour of leading order matrix elements, which in some respect can be 
understood as generalizations of the OPE. Identifying different classes of power-law corrections, 
they introduce a small number of universal parameters, e.g. effective values for integrals over 
powers of the strong coupling at low energies, which relate non-perturbative corrections for 
different event shape variables. 

The formalism has been worked out for mean values of some global event shape variables. 
From a measurement of the energy dependence of those mean values, a model independent 
determination of the strong coupling is obtained by fitting a s together with the non-perturbative 
power law corrections. The result of a first analysis of this kind is shown in Fig. 0. 



Figure 4: Measurement of a s 
based on mean values of global 
event shape variables measured 
between y/s = 14 GeV and y/s = 
133 GeV. 


Measurements of a s based on global event shape variables [|fi4| . |65| . [/]. (?/]. [G9|| from centre- 

of-mass energies between y/s = 10.53 GeV and y/s = 133 GeV are listed in Tab. Wherever 
available, the numbers are single experiment averages over several variables. Combining the 





















partial averages from energies in the range 10.53 GeV < yfs < 133 GeV into one global mean 
value finally yields a s (M z ) = 0.121 ± 0.005. 



Theory 

/GeV 

a s (ii) 

CLEO 

NLO 

10.53 

0.164 + 0.015 

TPC/27 

NLO+NLL 

29.0 

0.160 + 0.012 

MAC 

NLO 

34.0 

0.130 + 0.032 

MARKII 

NLO 

34.0 

0.153 + 0.032 

PLUTO 

NLO 

34.0 

0.108 + 0.038 

TASSO 

NLO 

34.0 

0.149 + 0.026 

MARK-J 

NLO 

34.0 

0.126 + 0.013 

CELLO 

NLO 

34.0 

0.144 + 0.026 

JADE 

NLO 

34.0 

0.162 + 0.043 

Average 

34.0 

0.134 + 0.019 

AMY 

NLLJET 

58.0 

0.130 + 0.006 

TOPAZ 

NLO+NLL 

58.0 

0.132 + 0.008 

TOPAZ 

NLO+NLL 

58.0 

0.139 + 0.008 

VENUS 

NLLJET 

58.0 

0.129 + 0.006 

Average 

58.0 

0.132 + 0.006 

ALEPH 

NLO+NLL 

91.2 

0.125 + 0.005 

DELPHI 

NLO+NLL 

91.2 

0.123 + 0.006 

L3 

NLO+NLL 

91.2 

0.124 + 0.007 

OPAL 

NLO+NLL 

91.2 

0.120 + 0.006 

SLD 

NLO+NLL 

91.2 

0.120 + 0.008 

Average 

91.2 

0.123 + 0.006 

ALEPH 

NLO+NLL 

133. 

0.119 + 0.008 

DELPHI 

NLO+NP 

133. 

0.119 + 0.009 

L3 

NLO+NLL 

133. 

0.107 + 0.010 

OPAL 

NLO+NLL 

133. 

0.110 + 0.008 

Average 

133. 

0.114 + 0.007 


Table 2: Measurements of the strong coupling constant from global event shape variables. The 
errors are the total errors of the individual measurements, which are dominated by the theoretical 
uncertainties. The averages were formed as described in the introduction. The results from AMY and 
VENUS are based on theoretical predictions calculated by the NLLJET Monte Carlo. The DELPHI- 
analysis at y/s = 133 GeV uses the second order prediction from perturbative QCD together with 
an analytical ansatz for non-perturbative effects. 


2.6 Jets from processes with variable Q 2 

One contribution to multijet production in ep-collisions is gluon radiation from a quark scattered 
off a virtual photon with large Q 2 . Quark and gluon emerge as two jets in addition to the jet 
from the proton remnant. The production rate R 2+1 of those (2+l)-jet final states is known 
to NLO for jets defined by the JADE algorithm. Tagging the scattered electron allows to 
control the Q 2 of the process and thus to establish the running of the strong coupling constant 
within a single experiment. First measurements have been published by the HI and ZEUS 
collaborations. [TOfl The conservative average of both results is a s (M z ) = 0.119 ± 0.013, which 
corresponds to ce s (fi) = 0.156±0.022 for a central scale fi = 19.6 GeV of the two measurements. 





















Measurement 

p/GeV 

Oi s O) 

Oi s {M z ) 

BjSR 

1.732 

0.320(|) 

0-118 ( 07 ) 

GLSR 

1.732 

0.260 (|Q 

o.no(8i) 

Rt 

1.777 

0.347 (37) 

0.122 (05) 

bb threshold 

3. 

0-217 ( 3 ®) 

0.110 (08) 

prompt 7 

4.0 

0.206(H) 

0.112 (l 2 ) 

DIS 

5.4 

0.200(H) 

0.116 (05) 

LGT 

8.2 

0.184(08) 

0.117(03) 

cc, bb decays 

9.7 

0.166(13) 

0.112 (06) 

Rnf 

18.0 

0.175 (23) 

0.128 (l 2 ) 

ep —Gets 

19.6 

0.156 (22) 

0.119(13) 

pp— > 66 +Jets 

20.0 

0.138 (f|) 

0.109(Q) 

e + e - fragm. 

36.0 

0.146(14) 

0.124(10) 

pp— >W+Jets 

80.2 

0.123 (25) 

0.121 (24) 

SM constraints 

91.2 

0.120(03) 

0.120 (03) 

event shapes 

< 133 


0.121 (05) 


Table 3: Compilation of measurements of the strong coupling constant. The total errors in the last 
two digits as given in parenthesis in most cases are dominated by the theoretical uncertainties. 


The error is dominated by the theoretical uncertainties. The published data still suffer from 
too low statistics for a convincing demonstration of the running of a si but improved results can 


be expected within the near future. |71 

Using radiative Z-decays, Z^qq'y, the running of the strong coupling can also be seen in 
e + e _ -annihilation by a single experiment. A first measurement with still rather low statistical 
precision has been presented by the L3-collaboration.[72| The strong coupling extracted at six 
different scales below the Z-mass was found to be consistent with the running as expected 
by QCD. The combined result is a s (M z ) = 0.119 ± 0.007, which corresponds to a s (p) = 
0.131 ± 0.008 for an average scale p = 50.8 GeV. 


3 A Global Average for a s 

A compilation of the results discussed in the preceeding section is given in Tab. [3]. The 
a^-values from the different types of measurements evolved to the scale of the Z-mass are 
displayed in Fig. Over an energy range from p = 1.73 GeV to p — 133 GeV all results 
appear to be consistent with one common mean for a s (M z ). The weighted average, assuming 
uncorrelated errors, is a s (M z ) = 0.1177 ± 0.0015 with % 2 /ndf = 6.7/14. The conservative 
average is a s (M z ) = 0.1177 ± 0.0037, where this error is biased towards large values through 
the inclusion of some less precise results. Restricting the average to measurements with error 
5a s (M z ) < 0.008, one finally obtains a s (M z ) = 0.1183 ± 0.0032. 

The value of the strong coupling constant at the physical scale of the respective measurement 
is shown in Fig. [|. The running in agreement with the QCD prediction is evident. Note that 
many different types of reactions with different numbers of active quark flavours are consistent, 
in agreement with the expectation that a s is flavour independent. 
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Figure 5: Measurements of the strong coupling constant from different sources evolved 
to the scale of the Z-mass. 


4 Summary 

In the past year significant progress has been made in the determination of the strong coupling 
constant by an improved understanding of power law corrections and perturbative higher orders 
on the theoretical side, and by precise new measurements from DIS, lattice gauge theories and 
standard model fits. All results are perfectly compatible with a common mean. Based on data 
covering the energy range 1 GeV < /r < 133 GeV, an average value a s (M z ) = 0.118 ± 0.003 
is obtained. The error takes into account the possibility of positive correlations between the 
measurements. From the RMS-scatter of the individual results an even smaller error would be 
derived. 
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Figure 6: Running of the strong coupling constant established by various types of measurements 
at different scales, compared to the QCD prediction for a s (M z ) = 0.118 ± 0.003. The open dots 
are results based on global event shape variables. 
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